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ABSTRACT 


The  state  variable  formulation  of  a  Foster  Wheeler 
ESD-Hboiler  is  developed  from  fundamental  principles. 
The  response  of  the  model  for  various  input  signals  is 
determined  using  CSMP  -l,the  IBM  simulation  language. 
The  sensitivity  of  the  model  to  various  coefficient 
values  is  noted  as  are  the  characteristics  of  various 
system  states  for  small  perturbation  values. 
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NOMENCLATURE 


As»Ay,A„  =  superheater,  riser  and  downcomer  cross-sectional 

2 

areas  respectively  (ft  ) 
C  =  evaporation  level  proportionality  constant 

(ft-sec/lb) 
Cs,Cy,CB  =  heat  capacitance  for  superheater,  riser  and  downcomer 

tubes  respectively  (BTU/lb-K) 
Ca  =  average  heat  capacitance  of  combustion  gases  times 

air- fuel  ratio  (BTU/lb-R) 
C;  =  neat  capacitance  for  f  eedwater  (BTU/lb-R) 
C^s=  average  heat  capacitance  of  superheated  steam 

(BTU/lb-R) 
Ss)  f*)  f  6  =  friction  coefficients  for  superheater,  riser  and 

2      5 

downcomer  tubes,  respectively  (sec  /ft  ) 

2 

g  =  acceleration  due  to  gravity  (ft/sec  ) 

h6  =  enthalpy  of  saturated  vapour  corresponding  to  P 

B 

(BTU/lb) 

"wb=  enthalpy  of  saturated  liguid  corresponding  to  P 

B 

(BTU/lb) 

KWB=   enthalpy   of    liquil    in   mud-drum    (BTU/lb) 

h-w  =   enthalpy    of    drum   and    downcomer   liquid    (BTU/lb) 

Hj3=   enthalpy   of    evaporation   corresponding    to   ?       (BTU/lb) 

B 

Kas=   air   cooling    coefficient   at   superheater   bank    (BTU/sec) 

Kttv,=  air   cooling    coefficient   at   riser   bank    (BTU/sec) 

Kgs'  Ks  =   heat-transfer    coefficients   from   combustion   gas   to 

superheater    tubes,    and   from   superheater   tubes 

o 

to   steam,    respectively    (BTU/lb-R) 
Kaa.Ky.)  Ki,    =    heat-transfer    coefficients    from   combustion    gas    to 


riser    tubes    and   from    riser  tubes   to    boiling 

4 

liquid   respectively     (BTU/lb-R)  ,  (BTU/lb-R    )  ,  (BTU/lb-R) 
Kft  =   evaporation    rate  constant    of   drum   liquid    (lb/sec-H) 
KT,KB=  constants  for    state    equations    of   saturated    steam 

(    R-ft  /lb)  ,     (ft  | 
LS)Lv>Lo   -   superheater,    riser   and   downcomer   tube    lengths 
respectively    (ft) 
M    =   mass    of   drum   liquid    (lb) 
Ms  =   mass   of   superheater   tubes    (lb) 
MB  =   mass    of   riser    tubes    (lb) 

2 

P&  =   drum    pressure    (lb/ft    ) 

2 

Ps  =  superheater  outlet  pressure  (lb/ft  ) 

2 

Pw  =   mud-drum    pressure    (lb/ft    ) 

&s =   heat-input   rate    from    tube   walls    into    the   superheated 

steam    (BTU/sec) 
0.35=   heat-input   rate   from    hot    gasses    into    superheater   tube 

walls    (BTU/sec) 
&b=   heat-input   rate   from    riser   tube    walls   into    boiling 

liquid    (BTU/sec) 
^3&=   heat-input   rate   from    hot    gasses   into   riser    tube   walls 

(BTU/sec) 
<i  =   feedwater    temperature    (E) 

e 

TB  =  saturation  temperature  corresponding  to  P   (R) 

B 

_  o 

Ts  =  superheater  outlet  temperature  (R) 

Tw  =  drum  and  downcomer  liquid  temperature  (R) 

Taw=  superheater  tube-wall  temperature  (R) 

Toy)=  riser  tube-wall  temperature  (B) 

T3S=  average  gas  temperature  at  superheater  banks  (R) 

t9b=  average  gas  temperature  at  riser  banks  (R) 

tc  =  comoustion  gas  temperature  entering  superheater  banks 

o 

(R) 

3 
VB=   volume   of   vapor   phase    in    drum    (ft   ) 
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Vy.  =   velocity   of    riser   mixture    (ft/sec) 
Vw  =   velocity  of    downcomer    water    (ft/sac) 

3 

V    =   total  drum   volume    (ft    ) 

Ws  =  steam  mass- flow   rate   at  the  superheater  outlet 

(lb/sec) 
Wf =   fuel    mass- flow   rate    (lb/sec) 
W;.  =   feedwater   mass- flow    rate    (lb/sec) 
Ww  =   downcomer   mass-flow    rate    (ln/sec) 
W   =   riser   mass- flow   rate    (lb/sec) 
Wtt=  air  mass-flow   rate    from   blower     (lb/sec) 
Wac=   chemically    correct    +50%  excess    air    rate    (lb/sec) 
We =   mass-evaporation   rate    from  drum   liquid   surface 

(lb/sec) 
WB =  steam    mass- flow    rate    from   drum   into    superheater 

(lb/sec) 
X   =   quality  of    mixture  leaving  riser 
*v  =   throttle   opening    (%) 
y  -   drum   liquid    level    (ft) 

3 

ph=  saturated    vapor   density   corresponding   to   P       (lb/ft   ) 

B 

3 

y°s  =  superheater  outlet  density  (lb/ft  ) 

3 

y°w  =   saturated   liquid  density    corresponding   to   P       (lb/ft   ) 

B 

/>    =   density   of    liquid   vapor   mixture   leaving   riser 

(lb/ft  ) 
Yj     =   evaporation    level    (ft) 
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INTRODUCTION 


Boilers  as  understood  by  marine  engineers  are  closed 
vessels  containing  water  which  by  the  application  of  heat  is 
converted  into  steam  at  any  designed  pressure.  This  steam 
is  then  used  for  the  production,  through  machinery,  of 
useful  work.  A  dynamic  model  of  steam  turbine  machinery 
consists  of  a  boiler  model  and  a  turbine  model.  Tne 
difficult  part  is  the  boiler,  as  a  load  change  causes 
variations  in  some  important  properties  such  as  toiler 
pressure,  temperature  and  drum  water  level. The  analysis  in 
this  paper  will  use  Chiens  dynamic  analysis  (1)  of  a  boiler 
as  a  reference. Chien  considered  a  naval  boiler  which  for 
purposes  of  analysis  was  divided  into  four  sections  namely  a 
superheater,  a  downcomer- riser  loop,  a  drum  and  a  gas   path. 

The  principles  of  thermodynamics, heat-transfer  and  fluid 
mechanics  were  used  to  describe  the  dynamic  behaviour 
corresponding  to  each  section  of  the  boiler  and  these  were 
derived  from  eguations  of  continuity,  energy, 
heat-transfer , and  momentum.  The  eguations  involve  partial 
differentials  as  well  as  nonlinear  terms.  These  eguations 
were  reduced  to  the  ordinary  linear  eguation  form  by 
applying  small  perturbation  and  difference  eguation 
technigues. Linear  eguations  thus  obtained  were  reduced  to 
ten  state  variable  eguations  and  solved  by  digital  computer 
technigues. 

Since  there  is  an  increasing  interest  in  toiler 
modelling,  tne  objective  of  this  thesis  was  to  develop  a 
comprehensive  boiler  simulation  model  in  a  form  useful  for 
modern  control  (i.e.,  mult ivariable  control)  analysis. 
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II.   BOILER  MODEL 


A.   BOILER  CONSIDERATIONS 


The  control  problems  of  high  pressure  boilers  have 
become  more  and  more  critical,  both  from  the  operational  and 
the  economical  points  of  view.  A  dynamic  analysis  is  the 
method  to  be  used  for  a  control- system  analysis.  It 
consists  of  a  complete  understanding  of  the  process  to  be 
controlled  and  the  effects  of  physical  and  chemical  changes. 
The  analysis  is  not  exact  by  any  means  but  the  results 
oDtained  should  be  in  good  qualitative  agreement  with  actual 
tests.  The  major  difficulty  in  boiler  analysis  is  the  fact 
that  the  whole  system  is  very  complex  and  contains  numerous 
variables. 

Chien(1)  considered  a  naval  boiler  which  for  purposes  of 
analysis  was  divided  into  four  distinct  sections  namely  a 
superheater,  a  downcomer-riser  loop,  a  drum  and  a  gas  path. 
The  detail  of  analysis  is  described  in  Part  B  of  this 
section.  The  basic  equations  used  in  the  analysis  are  those 
of  continuity, ener gy (heat-transfer)  ,  momentum  and  the  state 
equations.  These  equations  involve  partial  differentiation 
as  well  as  nonlinearities.  Generally, the  equations  have  the 
form 


•f  (x  >y,  2, )  =  o 


To  eliminate  the  nonLinearities   one   uses  perturbation 
theory,   which   effectively  approximates  the  response  of  the 
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system  to  small  signal  changes  about  a  chosen  operating 
condition.  Thus, the  equation  is  perturbed  about  its  steady 
state  operating  condition  to  give  the  linearized  form. 

Hence  it  can  oe  written  as 

h  ax  f  M  Ay  +  £f_Az+----  =  o 
ox     oy     dz 

The  perturbed  variables   are  AX,  AY,  AZ,   etc.   and  the 
partial   differentials   (  d£,  2>f ,  df,   etc.)   that   form  the 
coefficients  of  the  perturbed   variables   are   evaluated   at 
steady  state  operating  conditions.   This  technique  was  also 
followed  by  Whalley  (2)  in  modeling  the  same  boiler  plant. 
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SECTION  AX.      ELEVATION 


Figure    1   - 


A    CROSS    SECTIONAL    ARRANGEMENT    OF    THE    FOSTER 
WHEELER    D-TYPE     MARINE    BOILER 
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B.   BOILER  ANALYSIS 


The  boiler  studied  in  this  analysis  is  a  Foster  Wheeler 
D-type-  marine  boiler;  which  is  an  oil  fired, two-drum, natural 
circulation  unit  having  a  rated  output  of  28,800  lbs/hour  at 
350  lbs/in  gauge,  1200  F.  A  cross  sectional  arrangement  of 
the  unit  is  shown  in  Figure  I.The  following  assumptions 
apply  to  the  physical  simplifications  in  each  of  four 
sections. 


1 •   Superheater 

(a)  The  inertial  effects  of  superheated  steam  are 
neglected. 

(b)  The  superheater  tubes  are  assumed  to  be  a  single 
capacitance  with  restriction  on  the  drum  side  and  another 
restriction  on  the  load  side. 

(c)  Desuperheaters  are  not  considered. 

2 .   Down  comer  riser  loop 

(a)  Only  natural  circulation  exists. 

(b)  No  boiling  takes  place  in  the  downcomers. 

(c)  Vapor  and  liquid  velocities  in  the  riser  are 
identical. 

(d)  Heat-transfer  rates  to  the  boiling  liquid  from  the 
tube  wails  are  proportional  to  the  cune  of  temperature 
difference  between  the  wall  and  the  liquid. 

(e)  Steam  quality  is  uniform  in  the  riser. 

(f)  Liquid  temperature  is  always  the  same  as  the 
saturation  temperature  corresponding  to  drum  pressure. 

(g)  Downcomer  liquid  temperature  is  the  same  as  tne 
drum  liquid  temperature. 
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3.  Drum 

(a)  There  is  no  temperature  gradient  across  the  drum 
vapor  phasajand  the  temperature  is  always  the  saturation 
temperature   corresponding    to   drum    pressure. 

(b)  The  liguid  phase  has  no  temperature  gradient  other 
than   across   a    very   thin  boundary   layer    at    the   drum    surface. 

(c)  Evaporation  or  condensation  rate  in  the  drum  is 
proportional  to  the  difference  of  liguid  and  saturation 
temperatures. 

(d)  Feedwater    temperature   is  assumed    to   be   constant. 

(e)  Liguid-level  changes  due  to  bubble  formation  in  the 
drum   are   neglected. 

4.  Gas   Path 

(a)  The  air -fuel  ratio  is  assumed  to  be  constant. 

(b)  Temperature  of  combust  ion  gas  entering  superheater 
is  proportional  to  the  firing  rate. 

(c)  Waterwalls  are  lumped  with  the  riser-banks. 

(d)  The  heat-transfer  rate  at  each  tube  bank  is 
determined  by  the  tube  wall  temperature  and  the  average  gas 
temperature. 

(e)  Inertia  of  the  hot  gases  is  neglected,  that  is, 
velocity  changes  take  place  instantaneously. 

(f)  Delays  due  to  the  heat  capacitance  of  the  hot  gases 
are  neglected,  that  is,  temperature  changes  take  place 
instantaneously  in  combustion  gases. 

(g)  All  heat  transfer  is  due  to  turbulent  convection 
and  radiation. 

The  following  steps  are  taken  in  developing  the 
eguations  for  each  of  the  four  sections,  A  simple  schematic 
diagram  of  each  sub-section  is  included  at  the  beginning  of 
each  part  of  the  analysis.  A  brief  statement  of  the  physical 
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situation         is      included     under      the      headings.         The      dummy 

coefficients    such      as      a,       br      C,      D      ,      were      used      for 

1  1  11  11 

convenience^and   their    values    are   declared   in    the   Appendices. 

The    quantities  such   as    W    ,    h    ,    h    ,    etc.    are    the   steady   state 

B         S         B 

values   of   W  ,   h  ,   h  ,  respectively.   These  values  can  be 
B     S     B 

evaluated  from  the  original  unperturbed  nonlinear   equations 

by   setting  all  derivative  terms  to  zero  and  solving  for  the 
unknown  values  in  term  of  known  values. 
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Figure  2  -   SCHEMATIC  DIAGRAM  OF  SUPERHEATEI 
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5.      Superheater   continuity 

Mass      flow      in  -   Mass  flow   out   =   Hate   of   change   of   enclosed 
mass 

WB-WS         =    a.(LsAs/3s) 
The   perturbed  equation   is  : 

AWB-  AWS      =    LSASSA/S 
or, 

where      S      is      the '    Laplace   variable   denoting    differentiation 
with  respect   to  time 

o •      Superheater    energy    balance 

Heat   input    from   tube    wall    +    Heat   in    incoming   steam   -    Heat 
with   outgoing    steam   =    Rate    of   change    of  internal    energy 


Q+w.k-wh     =  i-fL sAapK    ) 
8         Bfc         3    s  rL  J 


The   perturbed   equation    is    ; 

AQs+WBAke+TiBAW-¥sAh3-hsAWs     =LsAs(/&SAKsfhsSAPs) 

Since,Ah&ao,  A^C^ATg     and  L_  As  S  A/s   =  AWB~AW§ 
After   substitution   and   rearrangement: 

(LsAsAs)siV(Sr^)iwB-Aas-¥^Ts  C«) 

or, 

a4SATs  +  asAWa  +  a6AQs    -    a?ATs 
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7.   Heat  conduction  across  superheater  walls 

Heat  input  to  wall  -  Heat  output  from  wall  =  Bate  of  change  of 
energy  stored  in  superheater  wall 

d 


Q   -  G   =I(HCT  ] 


9s    &    dt 

The   perturbed  equation  is    ; 

Aa9S-  Aa3  =m5cssatsw 

or, 

V^8  +  Vas=aioSATsw  (3) 


8«       Heat   input    to    superheater   wall 

The   following   empirical    formula   includes    the  cooling    effect 
if    air   supply    w^differs  from    ideal  air    supply   W      (whalley    (2)) 


ttC 


Via 
Wac 
The    perturbed  equation   is  '. 


^l.-".."*   (T,.-T„)-kai(l-^   )' 


Aa8s  '  lS  (T»»-1Vi)*,«ItS%f'(»V»T«)+2lUSs5?)". 


or, 


Vtt9«  "  ai2A^  +  Qx3(^9S-ATsw)+a4Awa  (4) 
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9«       Heat   transfer    across    superheater    wall 

For    turbulent    gas    flows: 

0.8 

o.     =   K  VI       f-!\.  -T  \ 
S  SB      I    8W      '  S ) 

The   perturbed   equation   is  : 

AS    -  |4.2Ks(TSw-Ts)AWa4.K3WB0-8(AT&w_ATs) 

ft 

or. 


ai5Aas   -  ai6AW*  +   ai7(ATsw-ATs) 


(5) 


10.   Heat  transfer  from  combustion  gas 

Average  gas  temperature  at  superheater  wall  is  T  where 

8s 


9S 


Tc- 


and   C   is  related  to  the  heat  capacitance  of  the  combustion 

3 
gas    by   ; 

Wo. 


c3  -c<^p 


where   C    is   constant 

c 

Thus  • 


a 


9S    -2C&(Wf+Wa)(Tc-Ts) 


The    perturbed   equation   is   ; 


Aft       =  2Cc[(w^Wa)6Tc  +  (fc-T96)A^  +  (fc-Tgs)AWa-(wf,Wa)ATgs" 
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aearranging ; 


or, 


2c  3*  '       <•'*.-'«.     3s' -    r     >-c     '5$'-"a,    n't    "a;u'9S 


ai8Aft9«-   ^AT<  +  a*oAVa2/\+VTqs  (6) 


1 1 .       Superheater   momentum 
Neglecting    the   inertia  of   the  steam  ; 

Since:  WB    =/>6AsV& 


^-sow* 


Substitute  :       X  f L 


i 


-s 


>     ~   9  V? 


The   perturbed   eguation   is    ; 

_2 
AP    -  AP  2<1»AWr  -i    J!"  Ap 

B       s        s  a  &A2  Jt 

Rearranging    ; 
or, 


a23A^+a,4AP3    +at.AA    =ae/PB  ^7) 


23 


*,/> 


T 


w">/>„ 


K>A 


WATER   DRUM 


Figure    3    -       SCHEMATIC    DIAGRAM    Of    RISER-DOWNCOMER    LOOP 
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12.   Riser  continuity  equation 

Mass  flow   into   riser  -  Mass  flow  out  =  Hate  of  change  of 
enclosed  mass 

The  perturbed  equation  is  : 


since 


1        „    JL     +     (i~x) 

y°         A        A 


then, 


V    -/(£-a)4x  +  #4^ 


If   the   starting   point    of  evaporation   in   the    riser   tube   is  Vj 
above    the    water  drum    exit   where   ; 


V)    -    1  -  cww 


then, 


A »  -  -CAW, 

Ay  ty  S  AP  can   be   neglected  since  >j  is   small 

Substituting   LP  and  AM   and    rearranging   ; 


o.3hYLYp  (4-  -  -4- 


)  SAX +  0.5  Ay.Ly.XKg    £     S^P; 


_2 


+  0.5  ( y5  +  ,ow)C  S  A  Ww 


=     AWW  -  AW 


or, 


«£7SAX  +  tf    8APfc  +  aa|SAW„    .  q30  aww  +  a3i  aw  (s) 
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1 3 .       Riser  momentum  equation 


Pressure   drop  -   frictional   losses   -    gravity    head  - 
exit   loss  -    monentum   flux   change  =  Inertia    force 


_■  p.\  - 


2. 
Wv. 


9  0°\)V         1CAAV): 


_        l_v> 


.W 


1    at 


Therefore  : 


9  2 

(p  -p  \    f  _l-/'iri:v+  3.  Y|jl  +      ww      _L  p  U    dw 


After   perturbation    and   substitution   for  Ap  as  done   before 
in   the   Riser   Continuity  Equation,    the    perturbed   equation   is: 


+  I|i  -L 


where  ". 


-2  "v 

ML 


I  A« 


F     = 

Bearranging  ; 


iih  + 1 


a/j 


Apw  -  apb  -  (V«,s)aw  +  bxAx-  -^a  Aww  +  -  ^ 


*  a;/w  w 


wnere 


b-  ■  "P£+«)/l4* 


a 


V      = 


26 


bv  = 


\A "  A) 


r  -z 


\a 


A1 


-   r  -2  / 1  , 


tl^+fJ+/V 


and   is   small    enough 
to    be  neglected 


Therefore  » 


APW-  APfi  =  ^AW  +  a^sAvUfc^AX-  *$*.   AWW 


or. 


a    APw  +  a     a  p.   „  a    A\ua    saw  +  q    *x  +  a37AW  (9) 
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14.       Hear   balance    across    riser   tube    wall 

Heat   from   hot   gas   to    riser    tubes-Heat    from    tubes    to    mixture 
:Rate   of    change  of    internal   energy   stored    in    tubes 


fi38-a6     -i(H»c.x„) 


The   perturbed   equation   is  ', 


Aft      ~    AQS 


M5C8S^BW 


or, 


a-  Aftq*  +    ft.  A  ft.  -  a,,  sat. 


38  V 


39  &  40  BW 


do) 
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1 5 .       Hi s er   heat    ba lance    equation 

Heat   input    to    risers-Heat    output    from    risers 
=Rate   of    change   of    internal    energy 

=  ^i[0.5A»(Lv-)))(/V,w)]+Av.(Lv^)(/^ja(h] 

From  the  continuity  aquation  ; 


3t 


w 


After  substitution  and  perturbation  ; 


k  =  K,,n  +  x  k 


since 

If    changes   in   latent    heat   ArU»are   neglected ,  then 


From   steam   tables    ; 

A  kn    =  Atw 
Therefore  : 

0.5  ^Lv(/+y5w)Kc8APa+  O.B  Av.LvC^+^J^     SAX  -  A<1B 

-(KWft-  K)a\-wkcapb-wTi^ax  +  vvwatw 

or, 

a4iSAP&+a42SAX+a4oAaB  =  a44AV'^a,3APB+a4eAX+a,7ATw  CiD 
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16.   Riser  heat  transfer  equation 

An   empirical   formula  that  accounts  for  turbulent  gas  flow 
is  adopted  and  gives  J 

The  perturoed  equation  is  ; 

Ad,,   =  "•^»(T..-fM.)  iWf  +  K3B5f  (AT„-4Ta„) 

>""(y'jH>4MT;VT:iTj 

Rearranging  ; 


ac 


0.6  0i4 


2Kfty(W&c-Wa)  3 

+  rS  ft      *    v  'ga  ATgh 


-   4Kv>TfcH  AT6., 


hence  : 


fll6K9i  Cjjia  -Taw)    a  w. 
+  -   0.4  t 

,      2Ka.v.(WaC  -via.)    AW 

—  £  a 

W 


or, 


a48Aa9&+  a49ATr      =  a5oATRW  fasiAwf +afi,Awa  (12) 
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17.   Riser  heat  transfer  to  boliinq  fluid 

The  heat  transfer  rate  is  assumed  to  be  proportional  to  the 
cube  of  temperature  difference  between  wall  and  mixture. 
Continuing  to  follow  Whalley (2) : 


a 


& 


VtBw-tj 


The  perturbed  equation  is  ; 


cr, 


Aa 


=  ay^w-TjCAT^-ATj 


a63Aa&=    as4ATBW    +a5SAT8 


(13) 


18.       Downcomer    momentum    equation 

Y,    Forces   due   to        pressure    diff erence-f rictional   loss 
+  gravitational   head-entrance   loss  =   Inertia    force 
Hence  : 

tj  b  ^  ij 


since  : 


^w    -/wAbVw 


2 
Ww 


(P  _P  )_  i^!w_  _ 

cr  p  D   A 2       <>'«    D      ooA/i 


.  £.  is. 


w„ 


K-p*)- 


fiEl->+o.5).J_lw2wf  pl  =  ii  i(ww) 

I      &A  A  ft..  A?  >    D        QAn   rlf 


3AA^ 


g*°  dt 
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The   perturbed  equation   is  ; 

APB  -  APW     -    (bw+Qws)AWw-L5A/3w 
where  : 


2Wy,      f   foLt 


aw  = 

Since   flow   is   incompressible    (  £>w  is   constant     both      in      time 
and   space), then     Apw  a0 
Therefore  • 


aseAPB  +  a87APvy    -  a58  Aww  +  Q59S  A  Ww  (14-) 

19.       Drum   st earn   mass   balance   equation 

Y.    Vapor    rates  to/from    drum    =   Rate   of    change   of      vapor      mass 

in   drum 

Thus  , 

Evaporation  rate  +  vapor  rate  from  risers-steam 

rate  from  drum  =  —  ( Va  Pa) 

at    J 
Therefore  : 

we  +  x'w-wB    m  i  (vBy>a) 

The    perturbed   equation   is  ; 

AWe  +  XAw  +  viAx-AwB  m,  V6s  Ay°&  +  /5&SAV& 
Substituting   for    ; 

AVa     ,   -AAy 

A /a    =>    KftAP6 

AWe    =    KefATw-ATB) 
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gives    ; 

VBKBSAPa-^aASAy  „  Ke(ATw-AT6)  +  XAW  +  WAX-AW, 


or, 


a6cSAP8  +  a6iSAy-a62(ATw-ATB)fa63Aw  +  a6+AX  +  aAvia  (is) 

20.   Drum  liquid  heat  balance  equation 

J2    Heat  to  and  from  drum  with  liquid  =  Rate  of  change  of  the 
internal  energy  of  the  liquid 
Hence  ; 

Internal  energy  of  riser  liquid  +  heat  transfered  in 
feedwater  -  heat  taken  in  by  downcomer  liquid  -  latent  heat 
released  with  evaporating  liquid  =  _d_(drum  liquid  internal 
energy)  «t 

Therefore  ; 

Cv(i-X)wTB+CiWiTi-CD^Tw-WeKBW  »i(CDMTj 

at 
But  Cr  ,  CL  j  Qb    -  i  • 

The   perturbed   equation   is   ; 

(t-X)f&AW  +  Cl-X)iiATft-V;TftAX  +  WlAT(i.Tuw{-W1,ATw 

-%AWw-WeAhBw-KMAWe      =       MSATW+TWSAM 

Substituting    for  ; 

AWe    =    KC(ATW-AT&) 


KcAP. 


AM     =   a/»w  A  y 
and    r.  »  constant, so     AT.   ^  o 
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Thus: 

i  Ci-X)TBAW4(i-X)WATB-WT8AX+TiAW.-WwATw-fvvAW,-WcKcAPa 

-KBWKe(ATw~  AT&)     =     MSATW  +  A/HfwSAy 
or, 

I  a66  aw  +  a67  atb  +  a6S  ax  +  a6g  awl *  a70ATw  +  an  aww  +  a72  apb 

I  +  a/3(ATW-ATB)  -      a74SATw  4-^_SAy  (16) 

21.       Dr u m   ligu id    mass    balance   ecjua tion 

£,    Liquid  input/output   rates   =   Rate   of    change   enclosed      mass 

of   liquid 

Hence  : 

feedwater    input   rate   +   riser    liquid   input    rate-downcomer 

liquid   rate-evaporation  rate   =   d  (  M ) 

Therefore   • 

W;  +   (i-X)  W-Ww   -  Wtf       =    i  (  M) 

The   perturbed   equation   is  ; 

AWj  +  (i-X)  AW-W  AX  -  AWW-  AWe     =   SAM 

Substituting   for  : 

AWe       =   Ke  (  ATW-  ATB) 

AM         -    A,PW  Ay 
gives: 

-Ke  ATB  +  A/wSAy  =    (l-X)AW-WAX-AWw-KeATw  +  AWi 

^6AT&+a77SAy   3   a84Aw+ag9AX+aaoAwM+a9iAVa9;Aw.    (17) 
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OV 


22.       Equation    from    throttle    valve 


\    '   M^^Tl) 


The   perturbed  equation  is    : 

AW     =  J&AXv  4--^  APS  +  |l4ATfi 

s  ox  H  ^Ts 

or, 


a9,    AWS      =    a94AXVf   a93APS  +  a9fe^ 


f  18) 


or. 


23.       State   equations 


AT.    =  1>  APa 


AT         :       a  AP 


-    KTAP5 


CIS) 


or, 


AA»trAP8     =kbap& 


iPi 


V«    =     a9*  A  Pe 


A  Tgs  .  ATg3 


C20) 


(21) 


ALw-^APa    . 


If. 


KCAPB 
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or, 


A  KBW      -  agq  a  p, 


99  «  rB 


(22) 


or, 


A  Ps 


-    -fz(Ts;A) 

=  Jis  AT   4-  lis  APe 


A  Pi        =    aio^Ts-+aioi  A/9 


C23) 


or, 


ATC        =    Kf  AWf 


ATC        -   &,„„  AWj 


102    ^ 


(24) 
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C.   STATE  VARIABLE  EQUATIONS 


The  following  procedure  is  an  algebraic  method  to  reduce 
the  twenty-four  governing  equations  to  ten  state  variable 
eguations#that  is,  in  standard  matrix  form  of  X(t)  =  AX(t)  + 
BU(t) 

Here  X (t)  is  the  matrix  of  system  state   variables   and 
U  (t)  is  the  matrix  of  system  inputs. 
From  (1) 


4A  *f,n  +  f;Aws 


or, 


AA  •  m*b  +  w*t 


(Li) 


From    (2) 


VTs 


-a  aw    -  a   aq  +  a  At. 


From    (5) 


(1-2) 


A  a     -    «L6Awa+.  £2A.T         jiaAT, 


(1.3) 


Substituting   Q   from  (1.3)  in  (1.2)  and  rearranging  ; 

s 


AT( 


a,6.a, 


i*±i7  +   Oz    AT  _    2»+  flu.ai6     >M  _   a<"  a^  At 
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or, 


AT     ~  b  At  -  b,  AUa  -  b_  At 


From    (3) 


(1.4) 


Oi.SATw  =  «,Aa3i+  a3  aqs 


From    (4) 


Ats,  =  ^iAa.ilAQ 
s*       alc      <js      ttlo       s 

A  a  .  =  ^  AWr  +  °i?  AT    .  ^3  at    +  a    A  w 
P       a-u       T       au      1s      au      5w         **      a 


From    (6) 


A  a     =  ^  AT  +ii?  AWi  +  iiiAwa+  £2?  AT 


Setting    (1.6)    =    (1.7)     and   rearranging*. 

[-] 


AT    -7 


flu       0^0 
0-U      0.1  J 


a  i3    02i 
ttii    ii* 


f  + 


0X3       02_2 
aU       CLlii 


_  Lillik'  Aw  +     ^ai3'    At 


an    aia 


£l_3        |£22 

an     u^ 


Substituting  AT  =  a  AWr    from    (24)    into    (1.8) 
rearranging  ; 


f*U    aio 


AT 


a 


Aia       Oj 

14         a  11      Ol8 


—i-  4.   "'"'" 


^13 

an 


3* 


&22 


aia 


AT  _ 


(1.5) 


(  i  •  6  J 


CI. 7) 


(1-8) 


and 


au" 


3 

J^L'  AW 


ittU    tti8 


Sw      /a         a    n  "    a 

"*13  ""22 


i^n     aia 
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or, 


AT.S    =    b6AWf+  t7ATsw-b3AWa 


(1.9) 


Substituting    ATa3from    (1.9)     into    (1.6)    and    rearranging   ; 
gs     lau        alt  )      t     \    alt       a*}     Srt     V  u        an    /      * 


or, 


AGl^.-   bgAWf  +  b10ATsw  +  baAwa 


(MO) 


Substituting      A  Gl      from      (1.10)       and       AQ..       from       (1.3)    into 
(1.5)    and  rearranging   ; 


0.1.0 


ATSw-|"-^-]^t  + 


'a8.bt0      ag.ai7 


\ 


al0  O-10'<^15j 


^+(^»a 


or, 


^M  ■  MWf  f  bi3  AVw  +  bi4  AWa+  b13  AWB  -  bi6  ATS 
from    (17) 


(Ml) 


A  y     ,  _  ^Zfi  ATR  +  ^15  AW  +.  ±ii  AX 


a 


77 


l77 


a 


77 


+.^AW,^^ii  ATWTi^  AW; 
a77  u.77  a77 


U.12) 
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Substituting     ATft    from    (19)    into    (1.12) 


8 


AY  ,  _  _!i_2l  apr   +  hi.  Aw+li.3  AX 

7  tt7;  5  0-77  »77 

4.  ^£  Aw,  +  ^2i  ATW    +  ^l  AW; 
a  77               a77  a  77 


or. 


Ay  =  -b17APB  +  bl8AW+  big  AX+-  b20AWw+b21ATw  +  b22AW:  (i.l3) 

From    (16) 

ATW  ,  -  1ZL  Ay  +  1"  AW  +  (!ill^)  ATfi  +  fiL  AX 


a 


7  + 


a  74. 


0-7+ 


a  7+ 


+ 


z:     i  "TtT"    w    a/*    w   a~ 


•74  a74 

« 


Substituting     Ay      from    (1-13)    and     ATR    from    (19)    into     (1.14) 
and   rearranging    j 


/\ 


.tw  _  ( *»■ bj7  |  Q97^67      gni2z3|  AP  +  |  *Si  -   a73,bia )  A\ 
a74     T      a 7+  a74  /  '    B   '    \a74  a7+    / 


,w 


,|^_  ^ZijAl    AX  +  [5a 
a™    /  a.74 


a78.b: 


a 


74 


74 


a 


74- 


Aww 


'S0+a73 


a 


a"-b«Wfc+'a69 


«».-bft,^ 


74 


a 


7+ 


a 


74 


a 


74 


AW; 


or, 


AT,   s   b 


*aAP»*  b24AW+^5^X  +  bi6AWw-^7ATwtb2yAWl     (i,l5) 


39 


From    (15) 


AP    __.?L'  Ay+^£  ATW_  iii£  Ax+l^AW+^i  Ax+^iiAW    (1-16) 
ft  a60     '      cuo  a60  a60  a6o  a6o      R 


Substituting     Ay     from    (1.13)    into    (1.16)    and   rearranging 


APD 


'«6i.bi7       Q-ti-tay' 


a  6o 


A  P.  + 


a. 


a 


60 


63 


e>        a6a 


i6i'bi8 


a 


AW  +    2>  - 


«64        a6i.bi3 


60 


l60 


/  a"  • bao  1  AW*   j. 


"2.6a         a6i-  bai  |    A_p 


a 


6.C 


a6c 


a  6o 


w 


flii.  b.j-i 


a 


6o 


»    AWi  +     -25    AW 


a 


6o 


or, 


AP&    -    b«  APB  +  Ko  AW  +  b31Ax  -b32^Ww 


+  b3aAT*  ~  b^^Wt  fb36  A 


wc 


Clbo 


AX 


Cl.17) 


From    (11) 

AX   .-^JAPj-  5+3  AaB  +iMAw,/-M  AP& 


a 


+2 


a« 


°-4z 


cua 


fi*J  AX  +  ^Z  AT, 

42  a42 


a 


(M8) 


From    (13) 


a  a     =  If*  AT_     +  i«L 


a 


BW 


53 


ATB 


a„- 


»3 


40 


Substituting    ATB     from    (19) 


ACL    .  f£i  At     +  fs-iisr  Ap 


aS3 


a 


5  3 


6 


C 1- 19) 


Substituting   A  Pp.  from   (1.17)   and   A  GL  from   (1.19)  into 
(1.18)  and  rearranging  ; 

a42    a42     a42.  aS3  '   B 


AX  -|;i"-  i+LSM   ^^'"^fa  f^L^l)^ 


a42 


\a4i         **2  /  I  a+2  a42  y  ^u       a42  y 


A' 


a42 


fa4i.b55>\ 


a 


42 


AwB  _ 


\ 


**-a«*|  AT, 

Q42.aS3 


&',V 


or, 


£  AX    .    b36A^-b37AW+b38AX+b33Aww 


+  b40  ATW  f  b4i  AW£  -  b42  AWB-  b43  ATBW 


From    (8) 


(1.20) 


AW   s-  °ZL  AX  _  i*?  AP,    +^3o   AW         Oaj  Aw 
a2g  a2g  azg        w      a29 


(1.21) 


Substituting   AX   from   (1.20)   and  A  P&  from 
(1.21)  and  rearranging  ; 


(1.  17)  into 


41 


AWV 


a*A*+^A9|Apa  + 


a 


29 


a 


29 


a27.b37       a28.b 


8-D3o      .     "31 


a 


+  ^£    AW 


29 


°*29  &29 


/ 


a19  U-23    i 


a29  U-23 


Aw, 


/ 


\ 


Q29  ai3      J       ^       1  a29  a  29       ■ 


a£9 

a27.b42     a2S.b33 


AW; 


a29  <*23 


AWR+. 


227-  l>48\     ^ 


/ 


a 


29 


BW 


or, 


AWW.  -baAPB  +  b43AW-b4,  AX  +  b47AWw-b48ATw 


+  b49  AWi+  bS0AW&  +  b5:  ATBW 


(1.22) 


Substituting     AWW  from    (1.22)    into    (14)    and    rearranging    : 


Pw  „f "»♦"»•  i>47  jAW       f°«,->a,,9.b44jA 


Ps  + 


I  &sa.OA! 


53    w43 


a 


57 


Aw 


a59.  b 


*L      AX  - 


fls«,  Oa* 


a 


57 


9  •  u4< 


0-67      . 


ATW   + 


al..t«  I  Aw 


a 


57 


a53-    bso 


a 


57 


AWbW^IAT. 


&Vi 
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or, 


APW   a   b32AWw-b53APfi  +  b34AW-b55AX-b56ATw 

♦  bS7Aw(+  bS8AW6+  b_9  ATBW  (1.2  3) 

From    (9) 

AW     -  ^i2  A  P    +   ^>  AP&  _  f3±  AW  _  £»i  AX  _  £sz  AWW 
"  Ss       W       a35  a3,  *33  Oa- 

Substituting     A  Pw    from    (1.23)    into    (1.24)     and   rearranging  ; 


a 


3S  a3.>5 


a»»         a,« 


3*      / 


a 


3  5 


a 


•Aw 


i5 


«3£_^53_  Q36 


\ 


a. 


a 


3.5 


AX  - 


fl32-t56 


a 


33       / 


*V 


^2.b57 


a 


AW, 


»32-b33 


AW. 


Vzz-bss 


35 


a. -is    y 


ATBW 


or, 


AW    ,   b6o  AWW  +  b61  APft  +  b62  Aw  -  b^  AX  -  bu  A7W 


'6a 


+  b65AW    +   b66  AW    +  b      AT 

1  B  0/  5W 


From    (10) 


(i-25) 


AT 


6W  - 


tt 


■40 


Aa 


^B         a4o 


a"    AQ 


(1.26) 
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From    (12) 

a4*A(^B     =    a5o^T8w+aslAW}  +  as2AWa-a43AT^  (1.27) 

Since     AT     =ATqs     from    (21), then    substituting     AT 
From    (1.9)    into    (1.27)    and    rearranging   ; 


Substituting     A  0LqFi  from      (1.28)       and      A&       from       (1.19)    into 
(1.26)    and   rearranging  ; 


♦l^Ur-*- MM^rN 


or. 


ATBW    =b68  ATftwt   b6gAWf  +   b7flAwa-b71ATsw  +  b72APB  (1.29) 

From    (7) 


AW     =    i£l  APft   _    £*£  AP    _   ^££  A/3  Cl<30) 


uu 


Substituting  A  Ps  from  (2  3)  and  A/>&  from  (20)  into  (1.30) 


Aw6  = 


<*23  &  23 


t..aice  \ 

b   i    a23  ;-■* 


V  «23     /     /* 


or, 


AW, 


=    b73APB-bHATs    -    b73   A/s 


(i-3i) 


From    (18) 


AWC   s   i£i  AXV  +    ^-s    APS    +  -5LSS.  ATS 


■95 


<X93 


l93 


(1-32) 


Substituting     Ap     from    (23)     into    (1.32) 


Aw.  .  ^1  axv  J 


a35.aloo  a96 


Cl93 


a93/      s  493       y 


or. 


Aws=  b76AXv  +  b77ATs+  b73Ay>s 

.Substituting     AW&   from    (1.31)       and      A  Ws       from       (1.33) 
(1.1)    and   rearranging    ; 

A/s  -  (  *>M\r  K'Ks)  A/t  ♦  (  b2  b77-  bt.b74)  AT£ 


(i.33) 


into 


^(W'b73  )APB  4  (b2.b7J  AXV 


or, 

Vs  -  Cu  A/,  +  C4-  ATS  +  di8  A  PB  +  DH  AXV  (A) 

Substituting     AW.  from    (1.31)    into    (1.4)    and   rearranging   ; 
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ATS    s   b4"  b75  AA  +  (  b3  +  b4-b74)  M&-  b3ATSw-b4b7.,APB 


or, 


ATS  =  d21  A/s  4  C22  At3  +  c23  AT3W  +  C2S  APB 


(B) 


Substituting  AWn  from  (1.31)  into  (1.11)  and  rearranging  '. 

O 

Atsw  =  -  biS.b75  Ay3s-(  b15.b74+bte)  ATS  +  b13  AT£ 


sw 


+  biS  .  b73  A  Pft  +  bi2  AWj  +  hu  A 


W, 


or. 


ATsw  -    tL31  A/>3   t  d32  ATS   +  C3a  ATSW  +  C3ftAPB 

+  D32AW^  +  D33AWa  ^C) 

Substituting     AWB  from    (1.31)    into    (1.20)     and   rearranging   * 

A*    =  b42:b73  A/s  +  K^7A  ATs+  k38AX-b37AW 
tb3g  AWW  -  b+3ATBW  +(b3b-  b4.2.b73)  AP& 
*■  b4-o  ATW  +  b4i  AWi 


or, 


AX    .  d4lA/>8    tC42  AT&  +  C^44AX  +  C45  AW  +  d46  AWW 


+  ^47  ATBW+Cl48AP6   +C49  ATW   +  D44AWi 


(D) 
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Substituting     AW     from    (1.31)    into    (1.25)    and  rearranging    ; 
AW  =  -b66-b73  AyOs_  b6€.  bM  ATs_  b63AX+  b62AW+  b60AWw 

+  b67  ATBw  +  (b66.b73+  b6i)  APA-b6+ATw+  b65  AW; 
or, 

0 

AW  =  e81A/>$  +  C52  ATS  +  ds4AX  +CiS8AW  +  dS6AWw 

^57ATBW+CJ5aAPft  +  d59  A7W  +DS4AW-  (E) 

Substituting     A  WR  from    (1.31)    into    (1.22)     and   rearranging  ; 

^ 

AWW  =,  -b50.b73  A/>s-  b5o.b74  ATS-  b46  AX+  b45AW 
I  +b47  AWW   +   b51AT&w+(bS0.b73-b4+)  Ap6 

-  b48  ATw  +-^49  AWt 


or, 


AWW  .  -C6i  A/>s  4  d62  ATS  +  d64  AX  +  d6s  AW  +  C66  Aww 

+  ^67  AT5W  +  cl68APB  +  C63  ATW+  DaAW;  (F) 

From    (1.29) 

0 

ATBw     =  -  b7i  ATSW  ♦  b^  ATbw   +b?2APft+bfc9AWf+b7oAWa 
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or, 

AT6W   =    d73  ATsw  +  ci77ATBW  +  ci78AP6H.D72AWf+D73AWa        (0) 
Substituting     AW&from    (1.31)    into    (1.17)     and   rearranging   \ 

A  i     -    -  Ka-  Ks  A/*  -    b35'  b74  ATS  4  bal  AX  +  b,0  AW 

\  ~b52   AWw+(b29  +  b35.b73)APa+b3:iATw-b^W. 

or, 

APB    =  C,tA/>6  +  C82  ATS  +d84  AX  +  dgsAw  +  C,&Aww 

+ci88  AP&  +  d83  ATw  +•  Dg4Aw£  (H) 

from    (1.  15) 

ATW  „    b^Ax  +   b£+ AW  +  b26AWw  +  b,3  APB+b27ATw^b23AW: 
or, 
I         ATW  -  C^AX  +dgsAW+c;gfeAww  +  ci98APB  +  dQ3ATw+D34AWi        (  I) 

From    (1.13) 

e 

Ay   =   bt9  AX+  bw  AW+bwAWH-bi7APB  +  baiATw4.b22  AWj 
or, 

|  Ay   .    C^AX  +  i10SAW  +  (i1HAWw-,Cl10?APfe+(i109AVDoAW,(j) 
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The  state  variable  equations  from  (A)  to  (J)  were  rewritten 
in  matrix  notation  as  shown  in  Figure  4.  The  state  equations 
were  solved  for  various  percent  step  changes  to  input 
variables  (  Ax  ,  Aw  ,  Aw  ),  using  the  IBM  simulation  language 
CSMP  -  IH  as  shown  in  the  next  section. 
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III.   COMPUTER  RESULTS 


From  Appendix  A  ,  it  is  seen  that  the  a  coefficients  are 
in  terms  of  steady  state  values.  These  values  were 
determined  from  the  data  in  Whalleys  thesis  (2)  which  are 
repeated  in  Appendix  E.  A  Fortran  IV  program  to  find  the 
a,b,C,D  coefficients  which  appear  in  the  CSMP  program  is 
shown  in  Appendix  F3and  the  calculated  values  are  shown  in 
Appendix  G.  Appendix  H  is  the  IBM  simulation  language  CSMP 
program  that  was  developed.  Only  input  variables  were 
changed  for  each  computer  run.  For  example,  results  for  a  5% 
step  change  of  Axv,  a  10%  step  change  of  AWjand  a  10%  step 
change  of  Aw:  were  calculated.  The  output  response  curves 
were  plotted  using  the  NPS  CALCOMP  PLOTTER. 
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A.   TRANSIENT  RESPONSES  FOLLOWING  A  STEP  CHANGE  IN  THROTTLE 
SETTING  OF  5% 


At  a  particular  steady  state  operating  condition  the 
throttle  valve  is  suddenly  opened  while  still  keeping  the 
previous  air  and  fuel  flow  rates  and  feedwater  flow  rate 
constant.  Because  of  bubble  formation  in  the  steam  drum 
liquid, the  steam  drum  water  level  (FIGURE  9)  swells  for 
about  40  seconds  and  then  shows  a  steady  falling  in  level  at 
higher  flow  rates  of  steam.  The  steam  drum  pressure  (FIGURE 
8)  immediately  begins  to  fall.  Following  these  sudden 
increases,  the  steam  mass-flow  rates  from  the  steam  drum  and 
the  superheater  (FIGURES  10  and  11)  show  slight  declines 
because  of  the  effect  of  the  decrease  in  steam  drum 
pressure.  The  superheater  outlet  pressure  (FIGURE  10)  acts 
in  a  similar  fashion  to  the  steam  drum  pressure,  that  is  ,it 
shows  a  steady  decline  to  a  new  steady  state  condition  but 
at  is  slightly  lower  than  the  steam  drum  pressure  because  of 
system  pressure  drops.  The  superheated  steam  temperature  and 
superheater  wall  temperature  (FIGURES  5  and  6)  both  decline, 
although  the  latter  temperature  is  slightly  lower.  The 
riser  tube  wall  temperature  (FIGURE  8)  and  steam  drum  liquid 
temperature  (FIGURE  9)  also  decrease  due  tc  the  throttle 
change.  The  superheated  steam  density  (FIGURE  5)  shows  a 
rising  characteristic  to  counteract  the  loss  in  superheater 
pressure  and  temperature  and  the  gain  in  steam  mass-flow 
rate  from  the  steam  drum  (FIGURE  11).  The  riser  and 
downcomer  flow  rates  (FI3URE  7)  are  reduced  to  a  lower  value 
but  steam  guality  (FIGURE  6)  is  increased  to  a  higher  value 
than  previously  to  preserve  the  energy  balance. 
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B.   TRANSIENT  RESPONSES  FOLLOWING  A  STEP  CHANGE  OF  10%  IN 
FUEL  FLOW  HATE. 


The  effects  due  to  a  change  in  the  fuel  flow  rate  are 
less  than  those  due  to  a  change  in  the  throttle  valve 
setting.  However, the  system  requires  a  longer  period  of  time 
to  reacn  a  new  steady  state  condition.  Because  the 
combustion  rate  is  increased,  the  evaporation  rate  is  also 
increased  making  both  the  steam  mass-flow  rate  at  the 
superheater  (FIGURE  17)  and  the  steam  mass-flow  rate  from 
the  steam  drum  (FIGURE  18)  increase  in  a  similar  manner. 
Initially,  density  in  the  risers  falls  quickly  in  response 
to  the  increased  firing  rate  which  causes  a  rise  in 
superheated  steam  density  (FIGURE  12)  .  An  increase  in  steam 
flow  from  the  steam  drum  without  an  increase  in  the 
feedwater  flow  rate  results  in  a  drop  in  steam  drum  watar 
level  (FIGURE  16). The  steam  drum  pressure  (FIGURE  15)  and 
superheater  outlet  pressure  (FIGURE  17)  rise  monotonically 
as  time  increases.  Also  superheated  steam  temperature 
(FIGURE  12)  ,  superheater  wall  temperature  (FIGURE  13)  , 
riser  tube  wall  temperature  (FIGURE  15)  and  steam  drum 
liquid  temperature  (FIGURE  16)  all  increase  with  increasing 
fuel  flow  rate.  The  abrupt  fall  in  riser  mixture  flow  rate 
(FIGURE  14)  is  probably  due  to  the  sudden  drop  in  the  level 
of  the  initial  density  change  in  the  riser. 
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C.   TRANSIENT  RESPONSES  FOLLOWING  A  STEP  CHANGE  OF  10%  IN 
FEEDWATER  FLOW-RATE 


The  main  effect  of  increasing  feadwater  flow  rate  while 
maintaining  throttle  setting  and  air-fuel  flow  rate  constant 
is  an  increase  in  the  steam  drum  water  level  (FIGURE23)  . 
Because  feedwater  is  added  into  the  steam  drum,  steam  drum 
pressure  (FIGURS22)  suffers  a  small  overshoot  followed  by  a 
decline.  As  a  result,  a  drop  in  both  steam  mass-flow  rate 
from  the  steam  drum  (FIGURE25)  and  steam  mass-flow  rate  at 
the  superheater  outlet  (FIGURE24)  occurs. Flow  around  the 
riser-downcomer  loop  (FI3URE21)  is  only  slightly  affected. 
Since  the  firing  rate  is  unaltered, the  riser  tube  wall 
temperature  (FIGURE  22)  and  the  steam  drum  liquid 
temperature  (FIGURE23)  decline.  The  superheated  steam 
temperature  (FIGU3E19)  and  the  superheater  wall  temperature 
(FIGURE20)  rise  because  of  a  drop  in  steam  mass-flow  rate 
from  the  steam  drum. The  quality  of  steam  (FIGURE20)  shows  a 
general  lowering.  The  superheated  steam  density  (FIGURE19) 
shows  a  decline  which  is  similar  in  form  to  the  drop  in  the 
superheater  outlet  pressure. 
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Figure     20    -      SUPERHEATER    WALL    TEMPERATURE    <     ATJW    )     VS    TIME 
AND    QUALITY    OF    MIXTURE    LEAVING    RISER    <      AX       )    VS    TIME 


70 


-O.DD 


-O.D0 


-D.  DQ 


52  m 


5 

C3 


-Q.DD 


G.tlO 


12.  DQ 


[X1D     J   *• 
16.00 


24.  DQ 


S.DD 


12- DD 


TIME 


cod 


12.  QQ 


16QQ 
CXia    J 


CXl3     J 
18- DQ 


24- DQ 


24.00 


E.aa 


L2.00 


TIME 


IS.  DO 

txin    1 


Z4.0D 


30.  DO 

a 
in 


1  o 

X 


30.  DO 


30-D33 




3 , 1 


<M     X 


3 


30.  DD 


Figure    21    -       RISES    MIXTURE    FLOW    RATE    <■      AW      )    VS    TIME    AND 
DOHHCOMER    LIQUID    FLOW    RATE    <       AWW    >    VS    TIME 


71 


-Q.QQ 


~a.ua 


-D.DQ 


6- DO 


12. 00 


CXia    J  * 

18-00 


E-OD 


TIME 


e-na 


12.00 


12.00 


IB.  00 

ixia    J 


[XiO     3 

18-  DO 


ixrn    J 


24.00 


24.00 


24.00 


24.  DD 


30.00 


;/"\ 

\ 

\ 

\ 

T 

3 


30.00 


3D.  DO 


3Q.QQ 
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IV.   CONCLUSIONS 


This  analysis  is  only  the  initial  step  toward  a  better 
understanding  of  a  naval  steam  generating  plant  both  from 
the  dynamic  and  from  the  control  points  of  view.  The 
analysis  concerns  a  boiler  system  which  is  divided  into  four 
parts.  The  major  difficulty  in  the  analysis  is  the  fact 
that  the  whole  system  is  very  complex  and  contains  numerous 
variables  which  are  extremely  unwieldy  to  manipulate.  Large 
numbers  of  relationships  are  non-linear  so  it  is  necessary 
to  generate  a  linear  form  by  the  methods  of  perturbation 
theory. After  setting  up  the  state  variable  equations,  the 
IBM  simulation  language  CSMP-EI  was  used  to  solve  for  the 
open  loop  transient  response  characteristics  of  the  state 
variables  due  to  variations  in  input  variables.  Throughout 
these  calcu  lations,  all  the  variables  appearing  in  the 
transient  response  curves  are  not  absolute  values. 
Rather, they  represent  a  small  incremental  change  from  some 
particular  steady  state  operating  condition.  Also, it  is  to 
be  noted  that  trends  in  some  of  these  results  differ 
markedly  in  form  from  those  appearing  in  Whalley  (2) .  Future 
studies  will  be  undertaken  to  determine  the  causes  of  these 
differences.  Meanwhile, the  benefit  one  may  derive  from  this 
basic  analysis  is  its  use  in  the  design  and  construction  of 
a  new  type  of  boiler  controller,  that  is,  a  multivariable 
control  system. 
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h 

fg 

= 

390.0 

h 
WO 

= 

0.0 

h 

= 

430.39 

K 
gB 

= 

3.6632 

K 

= 

7.5E-10 

r 

T 

a 

3000.0 

gB 


95 


bv 


ar 


A 
K 


BW 


W 
K 


=  930.0 

=  0.0 

=  0.059806 

=  890.0 

=  2.7014 

=  10.0 

=  0.3 

=  0.426 

=  40.0 

=  26.4 

=  -0.39157 

=  740.0 

=  894.38 

=  -1.7742 

■  405.0 

=  1680.0 
=  0.003 

=  500.0 
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APPENDIX    F 


FORTRAN    IV    PROGRAM 


$jc:b 


'T-'PLTC  i  T 

L S=45  .0 
A  S  =  0 . 2  3  9 

XbAL 

-**j- 

(A-Z) 

r(HOS=^  . 
CPS=0.5 
N  S  =  1  3  '>  0 


hB=1204 
WS=4.0 
M  S  -  L  9  0  0 

cs=o.?z 

KGS=0«g 

TGS=230 

-XSto  1.  2  3 

KAS=0.Q 


re- 


0  . 0 


WAC=7.0 

KS=4.68  4  2 
TS=1200 .0 

n  b  =  4  .  U 
COO. 37 

FS=6  J,.22 
RH0B=Q. 75a 

A  R  -  2  .  >  -> 
LR=1  ).0 
RHOlo. 5  2 

k  H  'Li  ,\  =  5  ?  .  r. 
X=  0.0  319  3.3 
K  h  =  n .  i )  n  c  o  ] 

C-C.J07637 

W  =  1  7  -j  .  •;•  7 

1  K  —  U  •  1  U  U  V  '-' 

DR=0. 146 
G=32.2 

MB=12 30  0.0 
CcJ=0.?.  ? 

K  C  =  0  .  0  G 
hFG=3  90 

HW 0=0.0 
H=430.3 

K G 3=3.6 


» Go 

.  n 


;<?  = 

TG& 

1  P  ■■ 

f.oh-\0 
=3  000. J 

--,  -;n  ,n 

K  -A  R 
KG  = 

=0.0 
D.05SE06 

FD  = 

L  U  = 

J  J  <  •  •  u 
2.7014 

I) . ..) 

L,H  = 

AD  = 

J.  3 

0.426 

4 ") « n 

A  =  2 
f-H= 

.6  •  + 

-0.39157 

7^0 »0 

r  iti = 3 


7   ij 
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■FBI 


M=1680.G 

kt=o.  in~. 

¥ F-5OO.0 

A I  =  1 .  0 
^2=-l .0 

A  3  = 
A  4-= 


A3  = 

A  6  = 


A3- 

Ac,= 


LS-A5 

I S*A  S*RHQS*CPS 

R  i-Hri 

-1.0 

-dS*CFS 


L0 

-1.0 

=  M  S  *  b  S 


A 1 1  =  1  .  0 
A12=(0.6*KGS)/(  Wl 


:.*)-m  rns-Ts**) 


A13=KGS*v>F**Q.6 


msfi 


K  A  S  /  K  A  e  ;: 


)•;-  H-.Jf-'-i&  ) 


A 1 5  =  1  .  0 

A16=(  J.  3JiiiS/WB**CUi  >  *<  ISVj-TS  ) 


i  /  =  <S~..-1**  J.  o 
A18=1.0/ (2.0*CG) 


4.4^ 


A- 


A20=TC-TGS 
A21=A2C 


«£  ^  -  -  (  s  F  t  ,■)  a  ) 

A23=2.0*FS*W37RHG8 
A_M=1  .0 


A25  =  -(  ['•S*..R**2/khL.j**- ) 
A26  =  l  .0 

,*s,i-r>*-  ?* 


^-7-=-0- 


.  K  a  L  i 


A<i3=0  .5*AR  *l.R*X*K3*  (  R! 
A29=0.5*C*{RH0  +  RFGi-»  ) 


(  t .O/XHHa-I  .J/Kh 


>/RHGP**2) 


A  -0=1 . a 

A31=-1.0 
A32=l.O 

A^j=-1.  J 

A34=2.0*«*(FR*IR/CR  +  L.5 )/(  G*AR**2*PHG) 

A3  6=(  1.0/RHH3-l.C/KhGVi  )  *(  (  W**2  /  (  G*AR**2  )  )*(FR* 
I ( RhO*#2*LR)  ) 

L  R  /  C  i- 

!  +  L 

.5)- 

A  J   1  -       £  •  U  T  /lift/    \  \1  I'M  f   '•'  ''"  £  T    '  r"  L  ,\  I 

A 38= 1.0 
A  3  9=- 1 .  0 

A  4  0 = m  a  -  C  B 

A'+1=0.5*AR*LR*KC*  (  r<r  J  +  R^n  -i ) 
A42=0 .  -">  *  <\  i<  * i  P  *  ri F  f-  *  (  F  i-  r  +  s  y-n ■•!  ) 

A43=-i.C 
A44=riWD-h 

A46=-W*FFG 

i  +  7-i.'i 

A  '-*  8  =  1.0 

A49=(  -KGS*WF#*0.6  )-<4.0*K.R*TG3**3) 

/\ -,o-|  - x  ."-  t) v  i»-F**0  t  *  )  -  (  4  . n*K  R * T P W **  ■* ) 

A 5 1=0.  6*KG3*  (TGb-TSiO  /V»F**0.4 
A5  2=2.0*KAR*{V,  AC-/.  A  )  /VAC**?. 

n  'J  J  =  I  .  U 

a  5  4=  3  .  0*K  ;*  (  Tt  'a-  T  l<,  J**2.0 
A5  5 =— A  5  '+ 

A  5  c  =  i.  .  0 

A57=-l  .0 

A^fi=?.c*i»n*(i^iU!  r/r:ntn.'i)/iin**7.n*r,*RHnui 

A5<5=!_0/ ( G*AO  ) 
A60=Vti*K  B 

AC2  =  K(- 
At  ;■  =  ,< 

-1.0 
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*  (  1  .  0  -  X 


AS7; 
A  6  8= 

A6rJ  = 


W*  ( 

-  W  * 


I  .  J-X  ) 


A7  J- 
A7M 

A  7  3  = 

_AL4i 


—  »•.  w 
-TW 


^nrz — 


—  M  (Z 
:-HB 


HCW*T 


*=4^* 


A9 
A9* 

-A-9- 


=  1.0 

=  o .  4 : 
o .  c  q  g  J  1 9 


A  96 
A97 

*.  ^  ■•■ 

A10 
AiJ 


0=  2  • 
1  =9 


000005 


]  C  J  .  J 


AH. 

Bl  = 

R  7- 


?=XT 

ai/a: 
a2/a; 


B3  = 
b4  = 

BJ- 

86= 
R?  = 


( A6#A17  )  / { A4*A 15 )  +  ( A7/ A4 ) 
( A  6* A16 ) / {  A4* A 1 b )  +  ( A  5/ A  4 ) 

(AC  »  A  1.7  i/U4*M-ttii- 


(  (  Al 
( A 1 3 


:,:  ft 


/  A 18  )  -  (  A  12/  Al  1  J  +  (  £20 /A  1 3  }  )  /  <  (  A 13/  Al  1 )-  (  A  22/A  1 8  )  ) 

I )  /  (  (  U  i  /  A  1 1 )  -  (  A  ?  2  /  A  1 1  )  ) 


(  A  I  4-  (  A  H  /  A  L  S  )  )  /  (  (  A  1 3  /  M  I  )  -  (  A  7  27  A  I  B  )  ) 

(  A 1 2  /  A  L I )  +  I  A  I  3  *  B  6  /  A 1 1  ) 

.=-U7*."  U/.U  i.  U-U  13/  All  ) 


Bb  = 
b9  = 
BIO 


Bll 

bL< 
-&*: 

Bl< 

B 1  b 


B17 


=Alt- 

= A  8  *  i 
-  (  A  6  ■■ 
=  A  3  * ! 
=  (A9  = 


=TAS? 
=  A  7  6  = 

-Ari-1  , 


(    1 

9/ 

11 
Al 


"AT 
A -3 
-A-Z- 


£*A  13 
A  1 0 

G  /  A  I  0 

/a; io 

6 )  /  (  A 


H7TA 
7/A7  7 

7 


/AH) 

)  ♦  (  (  A  j*,U7)  /  (  A  I  C* A  15  )  ) 


10*AJL5) 
lO*ALb) 


bl.9 
B20 


■5  ? 


A  B  9 
A  50 


A  7 
A7 


7 
7 

7 

1 

17// 

7TT- 


=A 

=  (  A7 

=  (  AG 
=  (  A7 


A7 


7T 
/A 
/A. 


4  1  +  (  A  9  ? 
A  /  5  ?  l 


Afc7/A7V)-('W7*A73/A74)HA72/A74) 
V  A  (  4  ) 


tiZ4 

fc25 


b27 

b23 


=  (  A7» 

=  (  A6' 
-(AC 
=  (  A  6 ' 
=  (  A  6  < 


+  A 
/A 

/A 
/A 


74* -( A7  5+B19/A74) 

74l-(  SJSS&ZOjLAJJl} 

74-( A75*B2l/A74 ) 
A75*B22/A74) 


7?)/ A' 

74)-( 

I  7/  A  t  C  ) — H'j^-ioy/A^C  ) 


TT77 
333 

r..,4 


6  0 )  -  ( 

'_21zl 

:J/AoJ 


A61*B13/A60 ) 

A  6  1  *  R  1 r  /  A  f.  3  ) 


=  (  AC' 


/A 


id)- 
2XAj 


(At>l*R21/A6  )) 


n  J  5 
B36 

p3  3 

r  3 


=A< 

=  (: 

=  (< 
=  (  1 


,5/ 
,45 

.4'j 


A.j 
/A 
-tr3r 
/A 
/  « 
7T 


+2  !-  (  A4  L*ti2 

3/A42 

i2)- (A4l*b3L/A4? ) 
'+  Z  )  -K  A  -.  L*;J2?/A42) 


7A4n  )  -  (  "i43-::AT.5*A9"/'\42/A  j3  ) 


X 


:  .7  -H  '  I 


B41 

I  i2 


k4l 
=AJLL 


A  4  2  *  b  2 

bi'J/A4; 
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B43 

344 


•■(  A4  3 
■■  (  A  2  7 


L'  4  5  -  t  ALT 
(  A  2  7 
( 


*A54 
*3  36 

■>  j  :  7 


)/ 

/A 


A  -'t  2  *  A  5  3  ) 

9)  +  ( A2  5*B29/A29) 

H — (  ,i2''**?'VA?-7  )  H  A3  I /A  20  ) 


346 

347 

-~~4~8~ 

e49 
BJLQ; 


?    5* 


•->.:>  . 


*B34 


=  (  A  2  8 
:  t  A  /  / 
=  A27*B43/ 
=  ( A5  3/A5  7 
■  t  A  '5  i-  /  A  9  7 

:A5  9* 

:  m  5  9  * 


/A' 
/A 

/a; 

ZA_ 


>)  K  A29*B31/A29  ) 

9)  +  (  A30/A29)-(  a  27*839/ A29) 

q  J  +  l  M'r  -*Hi3/A?9  ) 

9)-(A27*341  /A  29  ) 

si-(  A7.S*K-.'j/A2~») 


b51 

S52 

-3-§- 


A2 
)  + 


A59*B4"7/A5"7  ) 

5  9  ■:■-  l  t  'i  /  A  S  7  ) 


H54 
ti  5  5 : 


.4<-.  g 

B57 

ft  5  B 


=  A  3  9  ? 

=  A  5  9  * 


0<+5/ 

646/ 

C  +d  / 

349/ 

~-5_0Z 


A5  7 

A5  7 
TTTT 


A3. 


B6  2; 
B6  3= 


=A^9* 

=  {  A  3  2 
=  (A33 
=  (  A32 
■-(A3? 


B51/ 
"'5  3  2 

+M3 


<\5 

/A 


-(A3  7/A35) 


=  aj  r* 

:A32* 

■A12* 


*3  5  5 
EJ5T57 
3  5  7  / 


/A35)-( A34/A35 ) 
/A35 i  +( A36/A35 ) 


A-T5" 


do- 
B6 


a: 

A_3_ 


367  = 
668= 
B6  9- 


A  3  2  : 

(  A3 


359/ 
*  A  5  0 


f>»5U/4 


A3 
)/ 


A  40' 
4-4-fr= 


A  4 

-a-4- 


) +( A3 9*4 5 4) /( A40*A53) 

-)  -  {  A  1  S  »  A  4  9  * B  fc  )  /  (  A  A  Q*  HA  3  )- 


U  7  0  : 

671= 
-Pi 


(A3 

{  A3 

=  (A 


ci7  4 
B75 


=  (  A26 
=  A24 
=A24 


*A49 

*  A  5  1 


»/( 


A4  0*A48  )  +  (  43S*M9*BS  )/  (  A40*A4~  ) 
)/(A40*A43) 


T7: 

B77 

B-7-S 


'  A  9  4*7 
:  (  A  9  b 

:r\-5* 


:Ai.00 
;  A  1.  0  I 


/A  2 
/A2 


71/1 A40*A53  ) 

Set  -■-•  JS-'/A^j  ) 


A93 

*A100/A93 )+ ( A96/A; 

A-U1/A93 


CI  L 

LIZ 

-ebr 

C?l 

r  p  ? 


(S2* 
=  (32* 

:j4*d 

••■33+  ( 


tie) 

h!77) 


-(61 *  8 7  5 ) 
-( B1*B74) 


C2&; 
-C3JU 


7  5 

R4 : 


3/4  ) 


:  -  ,-j  r. 

:— B  4 

-  H  I  5 


^7: 


o_  2 

C2  2 
C3  5 
C4V 

C42 


=-(  31 

=  3 


5  *  B  7  4 }  -  B  L 


J  L  5  » 

:3i.P* 


F?75 
3  7  4 


L44--i3s 
C45=-&3  / 

C47=-S43 
C48=336-(84?.*B73  ) 

L  ■+  7-  1 1  +  J 

Co]  =-ta66  *B  75 

C52=-B66*3  74 

L34=-3c 3 
C55=862 

C5  7= jc  7 

C  5 8=  (  a  6  6  *  3  7  3 )  +  fi 6  1 

i   :■  ■,  -    lot 
C6l=- 65 0*375 

C  o  2-—  rj  '3  !J  -:: 3  74 

Co  5 

c  6  /. 


1TW 
=345 


C67  = 

Ct>b: 

Coa= 

r.73  = 


c51 

(  b  5  0 


* . ',  7  3  I  -  B  4 


— a  +  j 
-  3  7  I. 
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C  7  7  =  i3  6  9 

C  tl=oT7 
C8i=-b55*B75 


C84=B31 

Ce5=b30 
C36  =  - 


C  8  8  =  B 

Co  9-833 
L';4=S  >5 
C 95=82 4 

_CA_o=iL 


)  +  (  d  3  5  *  b  7  3  ) 


C98=623 
L99=8  2  7 

CI 94= R L  7 


C105=81 
C106=fi 20 

iu  .):;=-->  1.7 
CI  0 9  =  3  21 
D  U  =i> ' 


D3  2=B12 
033=814 


054=86 5 
K=B4  9 


DT2="3Ti9 
D73=B7  0 


10  0 


DO  4=  6; 
D104=L_ 

<RVTE — (  r  ,  U  0  i     44  j  J  ,  A  I  ")  1  ,  9 — 3>V'i, 


'  6  ,  g  7  7 ,  B  7  6 


FORMAT ( LHJL, 
19X.5H87 


!  X »  5  ri  A  i  0  0  = 
;  1 5  .  8  ,  5  X  ,  5 


,D  .  3  » 


sx 


hA101=,E15.8//, 


uj37^     =tF15.S,5X,5HB75    =,F15.8// 


:sx 


^m^MMIP 


■  mT1:^' 


LC42»C44,C4i>,C4fctC47.,C 

2 C 6 1 t  C  6 2 , C 04 

VM   l'nJ,;X,t»nL   LI 

=,F15.8,5X, 

=  .  E  I  5  .  3  ,  5  x 


l  :'.  i  i  C  2  2 1 5  g  5 ,  r.  >: 
ie,C49,C5X»C52tC54 


X*5HB78=,E15^E//| 


C3ttC^2f'CDj  >C 

,5  ,C56  ,C57»C58  ,C59, 


,C41, 


■u      TT 

19X,5HC21 


5K22 
5  bC3J 


T7TrTT5TTCTT 
=,F15.?,5X 


39X.5HC33 
^t9X  ,5hC42 

69Xt5HC49 


I15.fi  ,5X,5FC3r: 
lt>  .3,5X  ,5HC44 

•  i  •: .  r ,  3  x  ,  i  i  c  w 


tFl 


5TFTT:  :',  5 

5HC23    =»E  15.6//, 
X,5HC3?    =,F15.P//, 


XT5FCT3    =  .  E  E  5  .-5/7", 


=,E15.8,5X,5HC41 
,  E15.8,5X,5HC45    =  , 


79X 


,  j 


'C34 


P13 


jX,5HC51 

.  5X ,  5HC  55 


^,F15.F,5X,5HC5?     = 


o-X.VlLb  .' 
9S»X  ,  5HC61 
f,Ki  "L: U- 


=  ,E15  .3,  5X,  5I-C6 


,E15.S//, 

E  1.  5  .  8  /  /  , 
r  I  rj  ♦  C  /  /t- 

c  1 5  .  ?  /  /  , 
X",  5WC56    =  t  E  1  5  ♦  g  /  /  , 


,7,iV>)        C(: 


=  ,  t  L  5  .  cbX.b^L'j'.-    =TFT5TH7T, 
=,E15.8,5X,5HC64    =  ,C15.€/) 

;  r 'A ;,  ,  c  6  7  .  r  &  a  t  C  t c  t  r  7  -■ .  r  7  2  .  C  7  3  .  C  c.LJXil2-*£- 


>4,r£5_> 


lC86,C38,CB9,C94,C95,C96,C93,C9  9,Cia4,a%05,C106,C108,C109 
300       FORMAT  (If-    ,8X,5HC65    =  ,E15. 8 , 5X  ,  5HC66    *, E15.8 , 5X,5hC67    «-,El:>.8//, 

"13..-  ,oa  ,3'  &69    -  tri.^.Ct^X,bHC7.     ^Ll^.^/S- 


-ir;X<  5'T.  6  5 
29X  i5HC77 
3SX  ,5HC62 


V.X,  5*TC 


CX,5< 

a  .  3;-;r, 


J  S  X 

A.  c 


c 


=  » E 1 5 
=  fEl5 

=  i  E  J.  5 
=  ,  E  i  5 
=  .  E  I  5 


<»0C 


7SX, 3HCI 
:, -X  ,  5HC  I 

■■•■R'Tr — (- 


-T"f 


t1"   I 

B6-+- 


F0KMAT ( LHJ.SXt 

). r-  X,5HJ44    rrElj 

J  /  2    =TET5 


Z9X751 

3C.  X,  3. 


Lr 


, 5X,5HC7S    =,r1 .5.8,5 X.5HC 31 
,3X,3FC34    =,E15.3,5X,5HC85 


B",  ';X  ,  5HUTfc 
S,5X ,5HC95 

3  ,  n  X  .  5  b  f  1 9_ 


=,E15.8,5X,5HC96 


,E15.3//, 

= ,  E  .15.8//, 
=TET5TH77T 

,  F  1 5  .  S  /  /  , 


_EJL 


5_X  ,-^ri  ^=,MS.R//. 


,3FC I0ft= 


5X , 5HC108=, E  1 


//, 


=  //J 


i'A4,rc 


^r-* — ^m-'  4  ,  '.  -  -»  , 
H01 1    =,E15.3,5X ,5HD32    = 
6,5X,5H0  54    =tEl5.8,5Xf5H06^ 

T 


<t,  i?7?,  ")7?-,  On-i  r&94  ,9-W?4 

"X,5FC33    =,Fl5.t//, 


i/  ,  3 


=  ,<-  CT7Fr5T75HITM 


=  ,  E  1  5  .  ,r;  /  /  , 

=  ,  I    1  2  .  B77  , 


,3  <,5FPl« 


-  .  c 


1  3  .  S  /  /  ) 


END 
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APPENDIX    G 


FORTRAN    IV    PROGRAM    OUTPUT 


A  LOO=    D. 24000000b     )', 


/'1C1  =    0.9500COOOE 


B73    =   0.1415274CE-02 
B  Th    =    J .  +2  3UTJ  3  0  C t     7D 


b74     =    0.  3333n340F     00 


375    =    0.1 3 194440 E    33 

£72     =     C.16C49990E    01 


CI  I    =-G. L2436020E    J, 


C12    =-C.31416850E-01 


C  1 3 


13  15^21 OF-02 


C21    =     3.1713  3  75  JE     J  4 


C  2 


.16  0  = 


CI 


C35     =    3.5lc5fc340t:    01 


=~  L/  •  i  a  j  o  j>  j"  L  0 1  —  0 .1 
-^).5  4c65^!C3-0  I 


C31    =    D 


nJT^TiTE    u2 


C32    =    0.1  At 


5TCE  .5  J 
VI  OF  01 
Ob  =  -0.U94_1460F-02 
C4P    =    C.27G8C380E-05 


C33 


C  5  ■-: 


-0.4  7C.:-I1  c:F-C 


C44     =  --0U 


C42 
C46 
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